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I show how to compute the exact one-loop thermal correction to the free energy of a soliton. The 
method uses the effective potential as an auxiliary step to ensure that the soliton is quantized around 
the appropriate vacuum. The exact result is then computed using scattering theory techniques, and 
includes all orders in the derivative expansion. It can be efficiently combined with a calculation of 
the exact quantum correction to yield the full free energy to one loop. I demonstrate this technique 
with explicit computations in models, obtaining the free energy for a kink in 1 + 1 dimensions 
and a domain wall in 2 + 1 dimensions. 
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INTRODUCTION 

The effective potential is an essential tool in the analysis of the vacuum states of a quantum field theory. It provides a 
computational framework into which quantum Q| and thermal Q| corrections can be efficiently incorporated, allowing 
for analysis of phenomena including spontaneous symmetry breaking and restoration. A key characteristic of these 
problems is that they can be analyzed in terms of constant background fields. For background fields that are not 
constant, the effective potential is still useful as an approximate calculation, in which the contribution of the higher 
derivative terms in the effective action have been ignored. Once the length scale over which the field varies is no 
longer large compared to the Compton wavelength of the fluctuating fields, however, this approximation is no longer 
justified. In typical soliton problems, these two length scales are expected to be comparable. In this case, it is 
therefore desirable to replace the derivative expansion with an exact one-loop computation. 

For quantum corrections, such a procedure has been demonstrated in detail in |^, which implements ideas originating 
in [Q. In this case, it is essential to precisely implement the counterterms, which cancel the divergences of the quantum 
fluctuations. To obtain a meaningful finite result, these counterterms must enforce definite renormalization conditions 
fixed in perturbation theory. 

Thermal corrections are finite in the ultraviolet, since the contribution of a mode with energy to typically will fall 
like for large |w|. Indeed, they must be finite, since there are no additional counterterms available beyond 

those that have already been fixed at T = 0. While the absence of divergences means that the calculation of thermal 
corrections doesn’t involve the subtleties of renormalization, it introduces a different problem: Thermal corrections 
can modify the classical vacuum expectation value of the fields. While quantum corrections could in principle do so as 
well, there is a generally a counterterm available that can cancel this effect. Standard renormalization schemes specify 
that this counterterm be chosen so that the classical vacuum expectation value of the field is unchanged, which is 
represented diagrammatically by the vanishing of the tadpole graph. Such counterterms are T-independent, however, 
so they cannot cancel any further finite shifts in the vacuum caused by thermal corrections. 

Since analysis of the vacuum only requires constant fields, the thermal effective potential is well-equipped to analyze 
it. If a soliton exists at T = 0, then we expect that for small enough T, there should exist a similar soliton, built 
around the thermally modified vacuum state. The original T = 0 soliton configuration will now have infinitely higher 
free energy, since it approaches a value of the field that is not a minimum of the effective potential. 

Once we have constructed the thermally corrected soliton, the thermal effective potential gives an approximation 
to its free energy. But if we want the exact one-loop thermal free energy, accurate to all orders in the derivative 
expansion, we must instead sum the full contributions from all the small oscillation modes around the thermal soliton. 
At the same time, we must not double-count the approximation to this quantity already included by shifting the 
vacuum. We will see that the phase shift formalism of ^ provides an efficient mechanism for organizing this 
calculation. 
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This approach should be applicable for a wide class of held theory solitons, allowing for efficient numerical calcula¬ 
tion. Since we do not use any high-temperature expansions, we can continuously track the soliton starting from T = 0. 
To illustrate the calculation in a concrete way, we will choose a simple example, where much of the computation can 
be done analytically. We will consider the kink soliton in 1 -|- 1 dimensions, where the small oscillations potential 
is of the exactly solvable rehectionless Poschl-Teller form. We will also consider the corresponding domain wall in 
2-1-1 dimensions. Similar calculations have been considered in this model in Q. However, that work did not address 
questions related to the vacuum shift. 


EFFECTIVE POTENTIAL 

We will consider a model in 1 -|- 1 dimensions, with Lagrangian density 

where (f> is a real scalar held and m is its classical mass. Here C is a counterterm, which depends on the cutoh but 
not on T. At T = 0, the classical potential is 


= ( 2 ) 

and the model has a classical soliton solution 

TTIX 

4>kinii{x) = tanh — (3) 

which interpolates between the equivalent classical vacua (j) = ±1. Once we allow the temperature to be nonzero, we 
have to take care to ensure that we continue to consider solitons built around the correct vacua, which we will hnd 
using the quantum and thermal effective potentials. 

First, we consider the quantum effective potential, which we compute at T = 0. We use the standard method [Q to 
compute the effective potential and hx the counterterms. We choose a renormalization scheme in which we hold the 
location of the minimum of the effective potential fixed. This condition is equivalent to demanding that the tadpole 
graph vanish, and fixes C uniquely. Once we have defined the model in this way, the counterterm is fixed and cannot 
be changed when we consider T ^ 0. We also add an overall constant, independent of (j) and T, so that the value of 
the potential is zero at its minimum. This calculation yields the quantum effective potential 


VM = 


Stt 


C7(0) (^1-log^^ 


where the small oscillations potential is 




Following H], the thermal effective potential is given by 

rri poo _ 

= - log(l - 

TT Jo 


where T is the temperature. 

The total effective potential to one loop is then given by the sum 


H(<^) = W(</>) +fo, (</>) +. 


( 4 ) 

( 5 ) 

( 6 ) 


( 7 ) 


At T = 0, we have defined the theory so that the minima of this potential are at </> = ±1. But as T increases from 
zero, we cannot prevent the minima from shifting. For T large enough, the minimum will move to ^ = 0 and the 
symmetry will be restored. At that point, the soliton disappears. For T nonzero but well below this value, however, 
we will have a modified soliton. Our goal is to compute the difference between its free energy and the free energy of 
the trivial vacuum. 
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We can find the new minimum ipo{T) simply by minimizing eq. ( 0 ), which is easily done numerically. We must now 
build the soliton around this new vacuum. To accomplish this, we define a modified classical potential 


VM = w(<^) + - ^o(T)2) + BiT) 


( 8 ) 


where we have introduced artificial, finite counterterms A(T) and B{T). By choosing 


A{T) 

B{T) 


V^iMT)) 

200 (T) 




-W(0o(T)) = - —(0o(Tr 


1 ) 

If 


(9) 


we force Vc{(j)) to have its minimum at 4io{T), the same value of 0 as the full effective potential. Since we will only 
be interested in the difference between the free energy of the soliton configuration and the free energy of the trivial 
vacuum, we are free to add an overall constant independent of 0, even one that depends on T. We have used this 
freedom to force the value of the potential at its minimum to be zero. 

This modified classical potential is therefore identical to the zero-temperature potential, except that the artificial 
counterterm has shifted the mass of the small oscillations to 


ifi = \lm^ + ^^A(T) 

reflecting the vacuum shift caused by the thermal fluctuations. The classical potential is then 


Vc{4>) = 


2^2 


8A0o(r): 




( 10 ) 


( 11 ) 


and the kink solution is 


~ TTIX 

0kink(a;) = 0o(T) tanh — . (12) 

Since it is built around the correct vacuum, we will be able to use this soliton for calculations of the quantum and 
thermal corrections. At the end of the calculation, we will use its close relation to the original potential to subtract 
the artificial counterterms back out again, yielding the result in the original model. 


EXACT QUANTUM CORRECTIONS 

By considering the modified potential 14(0), we can now find the quantum corrections computed around the correct 
vacuum. The artificial counterterms A(T) and B{T) summarize the effect of the thermal contribution in shifting the 
vacuum. Later, when we compute the full thermal correction, we will subtract the contributions of the artificial 
counterterms and replace them by the full thermal corrections. The quantum calculation for 14(4) '^ow proceeds as 
a direct application of the techniques developed in ||^, |^ . We will outline this calculation, in the process introducing 
quantities that will be useful for the thermal calculation as well. 

We consider the small oscillations around the soliton, which are given by the solutions to the Schrodinger equation 

(-^ + t/(0(a;))^ 0fc(a:) = (fc^ -1 m^)0fe(a;) (13) 


where 


c/(0) = 4^"(<^) 

is the modified small oscillations potential. For the kink, we have 

t/(0kink(a:)) = (34kink(a^)^ - MTf) = ' 


(14) 


(15) 
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The quantum correction to the energy in the modihed theory is given formally by the sum in the shifts in the 
zero-point energies of the small oscillation modes, 




1 

2 


V 


+ Ect[ 4 >{.x)] 


(16) 


where the modes have energies Ej in the background and energies in the trivial vacuum. Ect is the contribution 
of the cutoff-dependent counterterm. To express this quantity more precisely, we will work in the continuum. Then 
the sum over modes is replaced by a sum over bound states plus an integral over scattering states, weighted by 
the difference in the density of states between the free and interacting cases. For a general soliton background with 
spherical symmetry, the spectrum can be decomposed in a partial wave representation into channels (. with degeneracy 
Df. In our one-dimensional example, these partial waves are just the symmetric and antisymmetric channels. In each 
partial wave, the difference in the density of states between the free and interacting systems is related to the scattering 
phase shift by 



Ptik) 

t-H I 

II 

1 

(17) 

and there are ni bound states, where 






(5^(0) = TTUi 

(18) 

which is modihed to 






(0) = 7r(n+ - i) 

(19) 


in the case of the symmetric channel in one dimension. In this case, there is a “half-bound” threshold state in the 
free background, if) = 1. As the name indicates, such states enter the sum over bound states suppressed by a factor 
of 2 1^. For the case of a charged particle, we must sum the phase shift for both signs of the energy corresponding to 
each value of the momentum k. 

Following 1^], we write the quantum correction as 

1 _, roo fjh fj _ rh~ 

4[(/)(a:)] =-E^j+y^ — +£ct[(t){x)\ (20) 


where u){k) = ^/k'^ + fh?, the ioj are the bound state energies, and S{k) = S-{k) + S+{k) = ^ logdetS(k) is the total 
phase shift in the potential of eq. dil)- The ^ term comes from the energy of the “half-bound” state in the free case. 
Since the phase shift falls like 1/k at large k, the integral in eq. (p0| ) is logarithmically divergent, and therefore should 
be considered as a function of a regulator, such as the dimension of spacetime |Q. The counterterm contribution also 
depends on the same regulator. Using Levinson’s theorem we have 

£q[(t){x)\ =+j —{u{k)-fh)—5{k)+£ct[(l){x)] (21) 

j ° 

and then we implement the counterterm contribution by expressing it in terms of the hrst Born approximation to the 
phase shift, 

^i(fc) = -^/ {U{(l){x)) - dx (22) 

giving 

£q\(j>{x)]^]^'^{u}j-m)+ ^{Cj{k) - m)^ (S{k) - Si{k)^ . (23) 

The Born term depends on the potential only through its integral over x, as it must since that is also how the 
counterterm depends on the potential. 
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In general the phase shifts and bound states can be efficiently computed numerically. Eq. 
exactly solvable Poschl-Teller reflectionless potential form. We have 

„ TO TO , 7 \ 3to 

d(k) = 2 arctan -p + 2 arctan — and oi(fc) = —— 

K 2/C rC 


where the bound states are 


(ho = 0 and uji = rh 


V3 


in addition to a “half-bound” threshold state at a; = m. Plugging these data into eq. gives 

fg[(?!'kink(a;)] = TO 

for the quantum correction to the modified potential. 


4v^ Stt 


however, is of the 

(24) 

(25) 

(26) 


EXACT THERMAL CORRECTIONS 

The one-loop thermal correction to the free energy in a background 4)q is given formally by a sum over the thermal 
occupation number of the small oscillation modes, 


8M{x)\ -^ 


-n\E^\IT _ 


Elog E' 


-n\E°\IT 


\n—Q 




\n—0 


= T ^ log(l - e-l^^l/^) - ^ log(l - . 


.-1^; 


l/T) 


(27) 


E° 


where we have taken the difference between the interacting and free backgrounds. The soliton has a zero mode, 
corresponding to a translation rather than an excitation. It should be omitted from this sum and handled by collective 
quantization, representing the motion of the soliton as a point particle. We will compute the classical energy plus the 
log of the full one-loop determinant with the zero mode omitted. This computation corresponds to an isolated soliton 
at rest. Ref. [D shows how to extend such computations to a dilute gas of solitons, which requires exponentiating our 
result and introducing additional factors of the classical action for each zero mode. 

As in the quantum case, we go to the continuum and represent the density of states in terms of the phase shift. 
Decomposing the scattering into channels i with degeneracy we have 


V- /v-'. A, A ^/t\ rdk^ /i_e—(D/T\ 

= T y: D, W log ( , _ ) - log (l - 0-”/-) + I - log ( ) 


J dk j 


— / dk 


(28) 


where we have used Levinson’s theorem in the second line.^ The prime indicates that the zero mode has been omitted, 
which also accounts for the mismatch with Levinson’s theorem resulting in the second term in the second line. In our 
one-dimensional example, this result becomes 


£t[()i(x)] = r j ^ log 


1 - 
I - 


dfc, /I - d(5(fc) 


i„g(i_o— 


dk 


(29) 


^ For fermions, where the background field is fixed externally, there is no vacuum shift. Thus we can compute the fermion free energy 
directly as s{ [(^{x)] = ^Ej + /o°° ^ log (l + This expression gives the thermal 

correction to the fermion free energy in a classical background. It is similar to the calculation used in for the thermal correction to 
the charge. 
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which is, as expected, a finite integral. 

Next we must modify this calculation to be able to consistently combine it with our manipulations of the T = 0 
classical and quantum contributions. There we introduced an artificial counterterm as a placeholder for the thermal 
shift in the vacuum. First we introduce these same terms into the thermal problem, so that we write eq. ( p^ in terms 
of the modified potential t/(0kink(a;)) built around the true vacuum (j)o{T). Next we must cancel out all the artificial 
terms we have added. We will accomplish this by imposing a consistency condition: We know that (f>o{T) is the 
local minimum of the full effective potential. By construction, it is a minimum of the modified classical potential, and 
remained a minimum when the modified quantum corrections were added, because we arranged to cancel the quantum 
tadpole. Thus if we cancel the modified thermal tadpole with a counterterm of the same form, we can be sure that we 
have compensated for all the previous manipulations. As with the quantum case, subtracting the contribution from 
the first Born expansion implements this prescription. It depends only on the background field through the quantity 
{U{(piainkix)) — m^), so it is indeed of the right form to implement this subtraction. Thus if we define 


£t[Hx)]=T 



/ 1 _ g 

V 1 - e-™/^ 




Jo ~ [ l-e--/T ) dk 


then the full free energy of the true kink system is given by 

d^[(l>{x)] = £c[(l>{x)] + iti^x)] + £q[(l>{x)] 



(30) 


(31) 


where we have defined the classical contribution 


£c[(l)ix)] 


Vc{4>{x))dx 


(32) 


which gives 

fh? 

£c[(t>kin]^{x)] = — (33) 

for the kink. We have thus made the unique choice such that the free energy is minimized at the true vacuum 
(j){x) = (j)o{T), and the value of the energy at this point is zero. 


DOMAIN WALL 


It is straightforward to extend this formalism to the case of a domain wall. As an example, we will consider the 
same cj)^ theory in 2 + 1 dimensions. This model was also considered in ©I- We now have the same kink solution, 
where the field is independent of the second spatial coordinate. Since there are no new divergences, we can continue 
to use the same renormalization scheme. In this case, the quantum effective potential is 


Vq{cj>) = 


127 r 


-U{(j))^U{(j)) + -mU{(j)) - -m^ 


and the thermal effective potential is given by 

l^t('/') = ^ f 'I'logCl - . 

Jo 

The method of [p^ gives the full quantum correction to the energy per unit transverse length as 


w 1 A- 2, r 2, 

= log 




ui{k)'^ d 
vn? dk 


[m-~s,{k)) 


(34) 


(35) 


(36) 


where the scattering data S{k), Si{k), and LUj are unchanged from the case of one space dimension, since the small 
oscillation wavefunctions can be separated into 


iljk,p{x,y) = V'fe(a;)e*™ 


( 37 ) 
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where p is the momentum in the transverse direction. Plugging in these data yields 

3to2 


IGtt 


(arccoth(2) — 2) 


(38) 


for the quantum correction to the domain wall. The analog of eq. (^) is now the thermal free energy per unit 
transverse length, 


27r 


' —OC 

TT 


y'log p"^ 


j - log (l - 


log ( ^ 


1 _ J dk 

where m{p) = \Jm? + ^ uj(^k,p) = + m{pY and uij(p) = . 


(39) 


EXAMPLES AND CONCLUSIONS 


We illustrate these results with some simple numerical calculations. Figure || shows the full effective potential in 
1 + 1 dimensions and 2 + 1 dimensions, with both plotted as a function of (f) at fixed T. 


V 




FIG. 1: Left: Real part of the full effective potential for the 1 + 1 dimensional model with T = 0.7m and X — 0.2m^, 

as a function of (p. Right: Real part of the full effective potential V{(f))lm^ for the 2 + 1 dimensional model with T — 2 Am and 
\ = 0.2m, as a function of </>. In both cases, the zero-temperature mass m sets the scale of units. Away from the minimum, 
the effective potential can become complex, which gives a cusp in the plot. 


Figure || shows the full free energy for the soliton and domain wall. In both cases, the plot shows the free energy 
of the extended object with the free energy of the trivial background subtracted from it, considered as a function 
of T, with the renormalized zero-temperature mass and coupling held fixed. The free energy cost increases as the 
temperature increases, as we would expect since the soliton and domain wall are low-entropy objects. 

This technique should be of general use in soliton problems at nonzero temperature. Since we have not relied 
on the high-temperature expansion, the calculation allows us to track the soliton continuously as we increase the 
temperature from zero all the way up to the point where the soliton melts away. In models with conserved charges 
|pd| , one could use this approach to compare the free energy of a soliton with a state carrying the same charge built 
on top of the trivial vacuum. It could also be used to study phase transitions and find percolation temperatures in 
statistical systems. One might also use the free energy to estimate rates of soliton formation after phase transitions. 
Finally, we have seen that this approach extends naturally to strings and domain walls using the techniques developed 
in ra. 
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FIG. 2: Left: Full free energy T[(l){x)\/m for the kink soliton in 1 + 1 dimensions as a function of T jm, with A = 0.2m^. Right: 
Full free energy per unit transverse length J^[(l){x)]/mL in 2 + 1 dimensions as a function ofT/m with A = 0.2m. In both cases, 
the zero-temperature mass m sets the scale of units. 
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